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Abstract 

In the present work it will be shown how the curvature of the outer cylinder 
affects Laser Doppler anemometry measurements inside a Taylor-Couette 



apparatus. The measurement position and the measured velocity are altered 



by curved surfaces. Conventional methods for curvature correction are not 



^ ' applicable to our setup, and it will be shown how a ray-tracer can be used 

to solve this complication. 

£>i! By using a ray-tracer the focal position can be calculated, and the veloc- 

i-£h ity can be corrected. The results of the ray-tracer are verified by measuring 

an a priori known velocity field, and after applying refractive corrections 

good agreement with theoretical predictions are found. The methods de- 

£>. | scribed in this paper are applied to measure the azimuthal velocity profiles 

in high Reynolds number Taylor-Couette flow for the case of outer cylinder 

^" . rotation. 

^sO ■ Keywords: Laser Doppler anemometry, Taylor-Couette, curvature, 

C*~) . ray-tracer, refraction, correction 
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1. Introduction 

K> ' A Taylor-Couette (TC) apparatus consists of two coaxial, differentially 

5— i ■ rotating, cylinders, see figure [TJ The annulus between the cylinders is filled 

with a working fluid; most commonly, as in our case, water is chosen. The 
apparatus has been used to study hydrodynamic instabilities, pattern forma- 
tion, turbulence, and was found to have a rich phase diagram with different 
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types of flow structures |ll4ll|. To get a deep understanding of these phe- 
nomena it is crucial to measure the local flow velocity. 




Figure 1: Left: Top view of TC apparatus, two concentric cylinders are rotating. Control 
parameters are the rotation rates an and ui , where the subscripts denote inner cylinder 
and outer cylinder, respectively. The inner cylinder has a radius of r^ = 20 cm, and 
the outer cylinder has an inner-radius of r — 28 cm. Right: Vertical cross section of 
a TC apparatus. The outer cylinder has an outer-radius of r e — 30.5 cm and a height 
L = 92.7 cm. The outer cylinder is made from optically transparent PMMA (Poly- (methyl 
methacrylate)) , and is attached to the top and bottom end plates. 



Measuring the velocity field inside a TC apparatus was done for a long 
time using intrusive measurements techniques, e.g. constant temperature 
anemometry 12H15I] and Pitot tubes [lg|. Though these techniques are ro- 
bust and proven to work, they are not ideal for measuring the velocity in TC 
flow. The aforementioned methods measure the magnitude of the velocity, 
not the individual components, and are directionally ambiguous using a sin- 
gle probe. Of course, one could use multiple probes 17[ to obtain the flow 
direction. Another problem is that they alter the flow under consideration. 
Though this is not an issue for non-recirculating setups, like an open-ended 
wind tunnel, it can be a severe issue in recirculating (closed) setups, e.g. 



a TC apparatus, a rotating drum, or a Rayleigh-Benard cell 18|]. For a 
larg e range of Reynolds numbers it is known that vortices will be shedded 
19( from these probes, either in form of a Karman vortex street or as a 
turbulent wake, depending on the geometry and Reynolds number. These 
vortices can survive a full revolution, which has been observed in rotating 
drum experiments [2CJ]. 

The TC setup used in the present work, the Twente Turbulent Taylor- 
Couette (T 3 C) [2ll423(]. distincts itself from other setups by many features: 
variable gap and radius ratio, precise temperature control, independently ro- 
tatable cylinders, and a fully optically accessible gap. The outer cylinder is 



constructed from 2.5 cm thick PMMA (Poly- (methyl methacrylate)) , which 
enables optical measurement techniques, e.g. Particle Tracking Velocime- 



try (PTV) |2J,|25j], Particle Imaging Velocimetry (PIV) |2Q,|27|, and Laser 
Doppler Anemometry (LDA) [281.1291]. These methods, by their very nature, 
will not disturb the flow under consideration. In addition these techniques 
are able to measure the velocity components and are directionally sensitive, 
such that they are capable of detecting flow reversals. The addition of seed 
particles is imperative for these techniques, and one should check if these 
particles accurately reflect the velocity of the flow, as discussed below. Addi- 
tionally, particles should not change the dynamics of the flow, in particular, 
some particles act as a surfactant in two-phase flows 3GH32I]. 



2. Laser Doppler anemometry 

LDA is based on the Doppler effect. The most common version of LDA, 
is a so-called dual beam heterodyne configuration [29J], see figures [2] and El 
In this configuration two beams are crossed and focused in the flow, creating 
an interference pattern. Seed particles, added to the flow, passing through 
the interference pattern will scatter light with a specific frequency. This 
light is then captured by a photo detector and converted to a current from 
which the Doppler shift can be calculated. Knowing the optical geometry of 
the setup one can directly calculate the velocity from the Doppler shift [291 ] : 

fd= T \Vk\ (1) 

where fa is the Doppler shift, A the wavelength of the laser, 6 the angle 
between the beams, and Vk the component of the velocity along ki — k2, 
where kj are the propagation vectors of the laser beams. To add directional 
sensitivity one has to frequency shift one of the beams, accomplished by a 
Bragg cell. More details about the B ragg cell, the fringe-model, and LDA 



in general, can be found in e.g. refs. [291. 1331]. 



Laser Doppler anemometry is a so-called absolute measurement method 
and therefore does not require calibration against a known flow. This, how- 
ever, does not mean that a measurement of velocity is error free. Any 
misalignment in the optical arrangement, and any imperfection in the lenses 
(e.g. astigmatism [3J]) will cause errors. In addition, any particle travel- 
ling through the beams prior to focussing can have adverse effects on the 
formation of a well-defined measurement volume. Similarly, any spatial in- 
homogeneity of the refractive index causes the focal point to shift, and the 



waists to mismatch in the measurement volume 35| . Furthermore, the num- 
ber of particles in the interference zone fluctuates; particles move in and out 
the measurement volume and induce noise in the collected signal. 





Figure 2: Left: The azimuthal and axial components of the velocity are measured, laser 
beams are in the green and blue planes respectively. Right: Vertical cross section showing 
two laser beams. The dashed lines are beams without refraction, the angle between the 
beams is denoted 8 a , where a stands for air. The solid lines are beams with refraction, 8 W 
is the angle between the beams in water, rt,n is the position of the laser head and tlh — / 
is position of the focus without refraction, while 77 is the real position of the focus. 



2.1. Curvature effects 

In most LDA applications the laser beams travel through flat surfaces, 
see figure El In this case, Eq. [T]can be simplified by invoking Snell's law: 



fd 



sin(6 w /2) sin(6j2) 



2\Vk\ 



A v 



A ; , 



where quantities with a subscripts denote quantities in air, and w in water. 
Equation [2] is only applicable if the interfaces are flat and the optical axis 
is perpendicular to those interfaces; it is only then that # a /2 is the angle of 
incidence and # w /2 the angle of refraction. The difference in refractive index 
is absorbed by the changing wavelength. So for the case of flat interfaces, 9 a 
can be obtained from the focal length and the beam separation, and together 
with A a , given by the laser, the velocity can be calculated from the Doppler 
shift (Eq. [2]). Note that the refractive indices of the container and water 
are irrelevant; they are not used in the calculation of the velocity. 

For the case of a curved surface, see Fig. El Eq. [2] does not hold. For this 
case Snell's law can not (easily) be applied in order to transform 9 W to 9 a . 
A prerequisite of calculating the correct velocity is therefore the knowledge 
of 6 W as a function of gap-position. 




Figure 3: Left: Typical geometry of LDA, equivalent to the vertical plane in the current 
application. The beams are passing through flat interfaces, and 9 does not vary with laser- 
head position. Right: Horizontal plane: laser beams are affected by the curved interfaces, 
and therefore 9 is a function of radial position. 



Most commonly the calculation of the velocity is implemented in the 
supplied software and implicitly assumes Eq. [2] to hold. For the case of 
curved interfaces this equation does not hold, and therefore the measured 
velocity has to be corrected by multiplying it with a correction factor Cq: 



Cg 



^</»,real 



n a sin{0 a /2) 



(3) 



^</>, measured Tl w Sin{U w / Z) 

where a subscripts denote quantities in air, and w in water, see also Fig. [2j 



3. Solutions 



roblem at hand is predominantly solved by mounting prisms (see 
to the outer cylinder of the TC apparatus, or by putting the 



371 ]) with flat windows. The 



The 
e.g. [3 

entire apparatus inside a liquid bath (see e.g. 
latter having two purposes: the liquid bath can act as a coolant and match 
the refractive index of the working fluid. In this way the beams travel 
through the outer cylinder with less deflection; this solution is, however, not 
perfect because of the finite thickness of the outer cylinder. Matching the 
refractive indices of the working fluid, the liquid bath, and the outer cylinder 
does solve the issue, but becomes cumbersome for large scale devices, or 
impossible if the studied fluid is a gas. The use of prisms is tantamount to 
the use of a liquid bath, and is also unable to fully correct for the problem. 
Furthermore, applying prisms is technically demanding once the outside is in 
motion. Theoretically one can derive the trajectories of the laser light. Ref. 
[381 ] derives these trajectories and even finds simplifications for the found 
solutions. This analysis is, however, not complete, as it only considers 1 



plane of refraction; the inside of the outer cylinder. In our application 
our cylinder is very thick compared to our measurement range; the gap of 
our apparatus is 8 cm versus a cylinder of thickness 2.5 cm. Here we will 
consider both interfaces, and calculate the trajectories by using a ray-tracer, 
and we will show that taking into account both interfaces is crucial for our 
experiments. Utilizing a ray-tracer has several advantages compared to a 
theoretical derivation for the present experimental setup (e.g. Ref. [38J): a 
theoretical derivation becomes cumbersome if one tries to find a formula after 
more than one refraction. The T 3 C system will be equipped with multiple 
outer cylinders to alter the gap-width, a ray-tracer is then more generalized 
and is able to handle multiple interfaces. The next section describes the use 
of a ray-tracer in order to account for the effects of the curved interfaces. 

4. Ray-tracer 

A 3D ray-tracer is built in order to calculate two parameters: the angle 
6 W for the green beam pair in the horizontal plane, and the position of the 
crossings of the blue and green beam pairs in the vertical and horizontal 
plane, respectively. For the blue beam pair (in the axial-radial plane), the 
curvature effects does not affect the flow velocity measurements, and there- 
fore Eq. [2] is applicable; the axial velocity does not require correction. The 
azimuthal velocity, however, does need correction. 

The ray-tracer is based on simple principles: starting at point p^ with 
direction kj, it checks which interfaces are hit for some t > at position 
Pi + ikj. The next point in the ray-trace can be defined from the interface 
that is hit first: Pi+i = pj + t m i n kj. The normal of this interface is calculated 
at position Pt+i, and is denoted s, where the hat means the vector has unit 
length. For the case of reflection the new direction is given by: 

k i+ i = k i -2(k i -s)s. (4) 

For the case of refraction, Snell's law: 

m(ki x s) = n i+ i(k i+ i x s), (5) 

is solved for ki + i under the constraint that it has unit length and in the 
plane spanned by s and ki. This can be implemented without the use of 
trigonometric functions, which can be troublesome in certain fringe cases, 
see the appendix for more details. 

Once the new position and new direction are found, the algorithm can 
be repeated until it exits the apparatus, or until it is absorbed by a surface. 



This algorithm has been applied to the geometry of our LDA and TC 
setup. Our focal length / = 0.5 m, our beam separation is 76 mm for the 
green beams, and 73 mm for the blue beams. The optical geometry of the 
TC apparatus used in the present work 2l|] can be characterized by 3 radii 
and 3 refractive indices, see Fig. Q] and [2] and Table [TJ 



Parameter 


Symbol 


Value 


Radius inner cylinder 


n 


0.20 


m 


Inner-radius outer cylinder 


r a 


0.28 


m 


Outer-radius outer cylinder 


r e 


0.305 


m 


Refractive index PMMA 


n PMMA 


1.49 


m 


Refractive index water 


71 

' ''water 


1.333 


m 



Table 1: Parameters describing the optical geometry of the presently used Taylor-Couette 
apparatus: the T 3 C [ll[. 



4 ■ 1 ■ Shift of focal position 

In this section the location of the foci {i.e. the measurement position) 
for both pairs of beams are calculated. The position of the laser-head tlh is 
varied, see Fig. [2j The focal position of the undisturbed beams is given by 
r LH — /• For each beam-pair the focal position rj is calculated as a function 
of rLH) see Figs. [2]and[H 

If the focus is outside the apparatus, i.e. tlh — / > r e , the focal position 
is given by rj = tlh — /• If the laser head is moved inward (decreasing 
r LH — /) the beams will first hit the outer cylinder at ry = r e = tlh — f ■ 
Moving the laser head further inward will cause the foci to lie inside the 
PMMA, and moving even further, inside the water. The blue beams are in 
the azimuthal-radial plane and refract differently from beams travelling in 
the axial-radial plane. The focus of the blue beams hits the inner cylinder 
(rj = rj) at tlh — / ~ 0.228 m, while the focus of the green beams hits the 
inner cylinder at tlh — / ~ 0.22 m. The distance between the two foci as a 
function of the radial position is depicted in Fig. [5j 

A pronounced shift of the focal positions is observed, see Fig. [5j the 
foci never coincide, and the maximum separation is 11.5% of the gap. The 
effects due to curved interfaces can therefore not be neglected. Furthermore, 
note that the foci do not even coincide at r = r ; this is due to the finite 
thickness of the outer cylinder, it is therefore necessary to consider both 
interfaces in the analysis. 
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Figure 4: The position of the two foci as a function of the position of the laser-head. 
Colors are in accordance with Fig. [5] The trajectory of the blue focus can be described 
by a piece- wise linear function, while for the green focus it deviates from a linear function 
due to the curvature of the interfaces. The foci diverge once tlh — / < r e . 
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Figure 5: The separation between the measurement positions is normalized with the gap 
width and plotted versus the focal position of the 'axial' beams. Near the inner cylinder 
the distance between the two foci is the highest and constitutes more than 11% of the 
gap. Due to the thickness of the outer cylinder there is non-zero distance between the two 
foci when focused at r — r . 



4-2. Beam angle correction 

For the case of flat surfaces (the axial-radial plane) the velocity calcu- 
lated by the supplier's software does not have to be corrected. Equation 



[21 however, does not hold in the azimuthal-radial plane due to the curved 
surfaces, and the velocity has to be corrected by multiplication with Cq, see 
Eq. [3J The refractive indices are known, and 9 a can be found from the focal 
length and the beam separation. 9 W can be found by calculating the angle 
between the focussing rays, see Figs. [2] and The correction factor can 
then be calculated, see Fig. [6l 



u 




r[m] 

Figure 6: Correction factor Ce (see eq. |3| as a function of the radial position r. The 
azimuthal velocity u^ has to be corrected: unreal = Ce(r) unmeasured- 



The velocity has to be corrected by 3% to 9%, depending on the radial 
position. Note that there is a finite correction (about 3% for the optical 
geometry of the current apparatus) at the water-PMMA interface at r = r , 
this is due to the non-zero thickness of the outer cylinder and taking into 
account both interfaces of the outer cylinder in the analysis. 



5. Experimental verification 

The results obtained from the ray-tracer can be verified by measuring a 
known flow state. The temperature of the TC system is kept at 20°C with 
water as the working fluid. Dantec polyamid seed particles (r scct j = 2.5 ixm) 
with a density of 1.03 • 10 3 kg/m 3 are used. One can estimate the mini- 
mum velocity difference Av = \v see( [ — vji U id\ between a particle v see d and 
its surrounding fluid Vfi u id needed for the drag force Fd rag = Girpr see( iAv 
to outweigh the centrifugal force F cent (r) = f 7rr seed 3 {p seed - p fluid) \- A 
typical velocity in the middle of the gap (r = 0.24 m) is v = 5 m/s, com- 
bining with the density and viscosity of water around 20° C, resulting in 
Av ~ 4 x 10 -6 m/s. This is several orders of magnitude smaller than the 



typical velocity fluctuation inside the TC-gap of order 1CP 1 m/s and hence 
centrifugal forces on the seeding particles are negligible. 

For TC flow a stable and well-known flow state is solid body rotation; 
the inner and outer cylinder are both rotated at a fixed speed to. After 
sufficient waiting the fluid will have a velocity u^ = cor, and u z = u r = 0. 
The experiment has been performed for three rotation rates oj (uj/2it = 1 Hz, 
2 Hz, and 4 Hz), where the azimuthal velocity has been measured at several 
radial positions and at mid height. Similar results were found for all three 
cases, Fig. [7] shows the results for the case of uj/(2it) = 2 Hz. 
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Figure 7: In red squares the uncorrected azimuthal velocity as a function of radial position 
is shown, and in blue dots the corrected azimuthal velocity. The black solid line is the 
theoretical flow profile u$ — cor, the deviation from this profile is plotted with green dots, 
the corresponding scale is on the right. The theoretical and corrected measured profiles 
are found to be in good accordance. 



The measured velocities are shown in red squares, after applying the 
beam angle correction the data points (blue circles) are found to agree with 
the theoretical flow profile within 0.75%. Any remaining deviation can be 
due to e.g. optical misalignment or imperfection, spatially inhomogeneous 
refractive index in the working fluid, or noise created by the amplification 
and digitalisation of the optical signal. 

6. Application 

Here the results are shown for the measurement of the azimuthal veloc- 
ity profile for the case that the inner cylinder is stationary and the outer 



10 



cylinder is rotating. This case has been studied before 12|, |l6(, and because 
the flow is laminar, even for high Reynolds number, any perturbation due to 
a measuring probe is likely to survive a full revolution. To accurately obtain 
the speed of the (undisturbed) flow, it has to be measured non-intrusively. 
Fig. [5] shows the results of three experiments having varying Reynolds num- 
ber (Re = uj r (r - Ti)/v). 
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Figure 8: For three rotation rates the azimuthal velocity is measured along the gap. This 
velocity is normalized by the driving velocity in order to collapse the data, and compared 
to the laminar velocity profile for infinite aspect ratio (dashed line). The profile is found 
to be nearly independent of the Reynolds number. 



For three rotation rates the obtained profiles are found to collapse over a 
decade of Reynolds number, but deviate from the laminar velocity profile for 
infinite aspect-ratio. Due to a finite aspect ratio of the setup the presence of 
the end plates will create a secondary flow, modifying the azimuthal velocity 
profile. The velocity close to the inner cylinder is currently unattainable by 
LDA due to reflections coming of the inner cylinder surface. These reflections 
create spurious, unreliable data and therefore the profiles are shown dashed 
in that region. The present system will, in the future, be equipped with a 
transparent inner cylinder, which will significantly reduce the reflections. 

7. Conclusion 

In order to measure correct velocities inside a TC apparatus using LDA, 
one has to correct for the effects of the curved interfaces. Not only do the 
positions of the measurement volumes depend non-trivially on the position 
of the laser head, but also the measured velocity has to be corrected. A 



11 



ray tracer has been used in order to calculate the position of the foci, but 
also to calculate the correction factor Cq as a function of radial position. 
We showed that for our application it is crucial to take into account both 
interfaces of the outer cylinder. The measurement positions do not coincide 
and the velocity has to be corrected even at r = r . Our ray-tracer is 
verified by measuring the velocity for the case of solid body rotation; good 
agreement with the theoretical prediction has been found. For pure outer 
cylinder rotation it is found that the velocity deviates from the laminar 
velocity profile for infinite aspect ratio due to the presence of the end plates. 
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Appendix 




t 



s 

Figure .9: Sketch of the rays in the plane of refraction. 

For each interface one has to solve Snell's law rii (lq xs) = 71^4.1 (kj+j x s) 
for ki + i. From the vector equation one can see that the solution is not 
unique; but the solution is unique if we require the solution to be in the plane 
spanned by s and kj (called the plane of refraction) . In this plane the solution 
can be rewritten as: njsin(#j) = rij+i sin(0j+i). ki + i can be decomposed in 



12 



to two parts; one in the direction of s and one that is tangential to the 
interface and in the plane of refraction, denoted t, see figure QJJ By simple 
geometry one can find that ki + i = sin(0j + i)t + cos(0j+i)(— is). The direction 
of t can be found by subtracting the normal component of the incident 
ray: t = lq — (Iq • s)s. From the definition of the cross product (|a x b| = 
|a||b| sin(a), where a is the angle between vectors a and b) and Snell's 
law, one can derive that sin(#j + i) = ^ k; xs , and using the Pythagorean 



identity one writes cos(#j + i) = \/l — sm(9 i+ i) 2 . Note that one does not 
need to calculate the angles Oi or #t+i> it is only necessary to know the sine 
and or cosine of the angles. The direction for the outgoing ray can now be 
found by substituting all the values: k; +1 = sin(#j + i)t + cos(#j + i)(— s). 
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